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ABSTRACT 

The transverse velocities of steadily rotating, non-axisymmetric patterns in 
flat galaxies may be determined by a purely kinematical method, using two 
dimensional maps of a tracer surface brightness and radial current density. 
The data-maps could be viewed as the zeroth and first velocity moments of 
the line-of-sight velocity distribution, which is the natural output of integral- 
field spectrographs. Our method is closely related to the Tremaine-Weinberg 
method of estimating pattern speeds of steadily rotating patterns, when the 
tracer surface brightness satisfies a source-free continuity equation. We prove 
that, under identical assumptions about the pattern, two dimensional maps 
may be used to recover not just one number (the pattern speed), but the full 
vector field of tracer flow in the disc plane. We illustrate the recovery process 
by applying it to simulated data, and test its robustness by including the 
effects of noise. 
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1 INTRODUCTION 



Over the past decade long-slit spectrographs have given way to integral-field spectrographs 
(IFS), which produce spectra oyer a fully sampled, two dimensional region of the sky (see e.g. 



Bacon et al. 



2001 



Thatte et al. 



2001 



Emsellem and Bland-Hawthorn 



20o3 ) . These spectra] 



* ssridhar@rri.res.in 

f sambhus@astro.unibas.ch 



2 Sridhar & Sambhus 

maps (also called the line-of-sight velocity distribution, hereafter LOSVD) contain impor- 
tant information on the flow patterns of non-axisymmetric features in galaxies and their nu- 
clei. It is widely believed that bars and spiral patterns in disc galaxies could influence galaxy 
evolution through their roles in the transport of mass and angular momentum. These pro- 
cesses are not understood completely, and IFS maps might be expected to play a key role i n 



the construction of dynamical models of evolving galaxies (|de Zeeuw 



2002: 



Emsellem 



2002). 



A limitation is that IFS maps provide information about radial, but not transverse, veloc- 
ities. It is not possible to recover the unmeasured transverse velocities without additional 



assumption s; a classic example is the model 



cular rings ((Rogstad. Lockhart. and Wright 



ingof the warped disc of M83, using tilted, cir- 



1974). However, the flows in non-axisymmetric 



features, such as bars, are expected to be highly non circular, and a different approach is 
needed. 



Tremaine and Weinber g (1984, hereafter TW84) considered steadily rotating patterns in 
flat galaxies, and showed how data from long-slit spectrographs may be used to estimate the 
pattern speed. Their method assumes that the disc of the galaxy is flat, has a well-defined 
pattern speed, and that the tracer component obeys a source-free continuity equation. The 
goal of this paper is to prove that, making identical assumptions about the pattern, IFS 
data can be used to determine not just one number (the pattern speed), but the transverse 
velocities, and hence the entire two dimensional vector field of the tracer flow. Like the 
TW method, one of the strengths of our method is that it is kinematic, and not based on 
any particular dynamical model. Our main result is equation © of § 2, which provides an 
explicit expression for the transverse component of the tracer current in the disc plane, in 
terms of its surface brightness and the radial current density maps on the sky. This formula 
is applied in § 3 to a model of the lopsided disc in the nucleus of the Andromeda galaxy 
(M31), where we also discuss the effects of noise on the data-maps. § 4 offers conclusions. 



2 THE RECOVERY METHOD 

An IFS dataset consists of a two dimensional map of the luminosity weighted distribution 
of radial velocities, the LOSVD. The LOSVD can be regarded as a function of the three 
variables, (X, Y, U), where X and Y are cartesian coordinates on the sky, and U is the radial 
velocity. The zeroth moment of the LOSVD over U is £ s k y (V, Y), the surface brightness 
distribution on the sky, and the first moment is F s ^ y (X, Y), the radial current density on 
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the sky 1 . Following TW84, we consider a thin disc that is confined to the z = plane, with 
x and y being cartesian coordinates in the disc plane. The disc is inclined at angle % to the 
sky plane ii = 0° is face-on, and i = 90° is edge-on), with line of nodes coincident with the 
x-axis. It is clear that the sky coordinates, (X,Y), may be oriented such that the X-axis 
and X cLXIS £1X6 coincident. Then (X,Y) = (x,ycosi) . 

The non-axisymmetric pattern of the tracer is assumed to rotate steadily at angular rate, 
Q p z . In this frame the continuity equation for the tracer brightness assumes its simplest 
form. Let r be the position vector in the rotating frame, E(r) the tracer surface brightness, 
and v(r) the streaming velocity field in the inertial frame. An observer in the rotating frame 
sees the tracer move with velocity, [v(r) — Q p (zxr)] . If the tracer brightness is conserved, 
E and Ei> must obey the continuity equation, V • [E (v — Q p zxr)] = 0. Cartesian 
coordinates in the rotating frame may be chosen such that they coincide instantaneously 
with the (x, y) axis; thus r = (x,y) andf(r) = (v x , v y ) . In component form, the continuity 
equation reads, 

ox ay \ oy ox J 

which is equivalent to equations (2) and (3) of TW84. The quantities, E(x, y) and E(x, y)v y (x, y), 
can be related directly to the observed surface brightness and radial current density maps: 

E(x,y) = cosz E sky (X, Y) , (2) 

E(x,y)v y (x,y) = cot i F sky (X, Y) . (3) 

Henceforth E(x, y) and E(x, y)v y (x, y) will be considered as known quantities. The unknowns 
in equation (P) are Vt p and E(x, y)v x (x, y). Below we prove that both quantities may be ob- 
tained by integrating over x . We will assume that E(x, y) , E(x, y)v y (x, y), and (the unknown 
quantity) E(x, y)v x (x, y), all decrease sufficiently rapidly with distance, such that all the in- 
tegrals encountered below are finite. 

Integrating equation (PJ) over x from — oo to x, we obtain, 
d f x 

E(x,y)v x (x,y) = -— J ^ dx (Ev y - fi p x'E ) {x , y) - Q p yE(x,y) , (4) 

where we have used E(— oo.y) = 0, and E(— oo, y)v x (— oo, y) = 0. We must also require 
that E(+oo, y) = , and E(+oo, y)v x (+oo, y) = . This leads to the condition, 



1 The mean radial velocity is then given by U(X, Y) = (Fgky/Ssky) : the contour map of U(X,Y) is often refered to as a 
"spider diagram" . 
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— J dx (Ei^ — VtpxT?) = . (5) 

Since the integral in equation (J5J) is independent of y, we can infer its value at large values 
of \y\. Therefore the integral itself must vanish, i.e. 

/+oo r+oo 
dxxT,(x,y) = / dx E(x, y) v y (x, y) , (6) 
-QO J — CO 

for any value of y. This will be recognised as the key relation that TW84 employ to determine 
the pattern speed (see eqn. 5 of their paper). As is clear from our derivation, the real 
significance of equation (0) is an eigenvalue of equation (jlj). In other words, it provides a 
consistency condition that E(x, y) and E(x, y)v y (x, y) must satisfy, if E(x, y)v x (x, y) is to be 
given by equation (jlj). Using equation @, we can rewrite equations (0) and (JH), such that 
fl p and E(x, y)v x (x, y) are expressed directly in terms of observed quantities 2 : 

/+oo r+oo 
rfXIS sky (X,y) = / dXF sky (X,Y), (7) 
-oo J — CO 

(Eu,)^) = -(cos*) 2 ^ " ^ x ' Ssk y) (x , y) " ^s sky (x,y). (8) 

In the next section, equations (JTj) and (jHJ) will be used on the simulated data of Figure (JTJ), to 
enable recovery of the entire two dimensional flow vector field of a steadily rotating, lopsided 
pattern. 



3 APPLICATION TO SIMULATED DATA 



We used as data, simulated observations of a numerical model of the stellar disc in the nucleus 
of the Andro meda galaxy (M31). A brief ac count of the model is given below, and the reader 



is referred to 



Sambhus and Sridhar! ()2002l hereafter SS02) f or details. The nucleus of M3 1 



is believed to harbour a Super Massive Black Hole (SMBH) f Ko rmendv and Bender 



surrounded 



jy a d ense stellar disc which appears as a lopsided, double-peaked structure 



(LaueretaL 



1993). The two peaks are sepa rated by about Of! 5 , with the fainter peak almost 

re dynamical cen tre of the 



coincident with the location of the SMBH ( LaueretaL 



1998J). T 



Tremaind (|l995[ ) proposed 



galaxy lies inbetween the two peaks, about Of!l from the SMBH. 
that the SMBH was surrounded by an eccentric disc of stars, whose orbital apoapsides were 
aligned in a manner that gave rise to the lopsided peak in the density of stars. Our input 



2 In principle, the value of Q p given by equation J3 should be independent of Y. However TW84 recommends multiplying 
equation by some weight function, h(Y) and integrating over Y , to obtain an estimate of O p as the ratio of two integrals — see 
equation (7) of TW84. 
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Figure 1. Simulated data from a model disc inclined at 51?54: two dimensional distribution of (a) Surface brightness (Sgw) 
and (b) Radial current density (-F s ky) of the model disc. The images have been smoothed with a circular Gaussian beam 
(cr = O'.'l). The contour levels are arbitrary, but separated uniformly in the values of T, s k y and F s k y , respectively. In (b) the 
black and white shadings correspond to negative and positive radial current densities, and the dashed line is the zero radial 
current density curve. In both maps, the line-of-nodes is along the x axis. The axes scales are in ("). ("Data" taken from SS02) 



model is a dynamical m odel of this eccentr ic disc that was constructed by SS02, based on 



the HST photometry of 



Lauer et al 



( 19981 ). The model consisted of about 230,000 points 
distributed on a plane. Each point ("star") possessed five attributes: luminosity (or mass), 
location in the plane, and two components of velocity. The lopsided pattern formed by 
these points rotated steadily about an axis normal to the plane with a (prograde) pattern 
speed equal to 16kms _1 pc _1 ; thus the model satisfied all the assumptions used in § 2. SS02 
estimated an inclination angle i = 51?54, and we use this value while projecting the model 
disc to the sky-plane. To obtain a smooth distribution, we "observed" the model with a 
circular Gaussian beam of a = O'.'l . Figure (0) shows the surface brightness (S s k y ) and the 
radial current density (-F s k y )- The line of nodes is along the x-axis. 

We computed the integrals in equ ation Q, using E s k v and F^ v from Fi gure ((H), for 



eleven different values of Y. Following 



Gerssen. Kuiiken. and Merrifield (1999), we plotted 



the eleven different values of one integral against the eleven values of the other integral. The 
slope of the "best-fit" straight line (in the least-square sense) gave an estimate of fl p sin i. 
Using i = 51?54, we found that Q p = 15.11 ± 0.47 km s -1 pc _1 . We used this value of fl p to 
compute the right side of equation (JHJ). After deprojection using (X,Y) = (x,ycosi) we 
obtained the x-current density, (St> x ), the isocontours of which are displayed in Figure (2a) 
as the continuous curves. For comparison, we also plot the isocontours of (Ev x ) from the 
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Figure 2. Isocontours of the i-current densities (a), and j/-current densities (b) in the disc plane. The contour levels are 
equally spaced in current units. The continuous and dashed curves correspond to the recovered and model current densities, 
respectively. The black and white shadings correspond to negative and positive current densities. The velocity field of the input 
model is shown in (c). The recovered velocity field is close to the input model, and we show only the residuals (i.e. recovered 
minus model velocities) in (d). All axes arc in ("). 

input model in the same figure as the dashed curves. In Figure (2b) similar plots of (T,v y ) 
are displayed to make the point that, in practice deprojection can also give rise to errors. 
It is traditional and useful to look at the velocity field, instead of the current density field. 
The velocity field is obtained by dividing the current density field by and we may 

expect this process of division to give rise to errors, especially in the outer parts where 
£ is small. To quantify the errors, we computed the residual-map, which was defined as 
the difference between the recovered and input x-velocity maps. The S-weighted mean (R) 
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Noise 


Hp (error) 


R (error) 


ctr (error) 


1% 


15.29 (0.79) 


1.14 x 10~ 3 (1.75 x 10~ 4 ) 


8.02 x 10~ 2 (1.54 x 10" 3 ) 


5% 


17.67 (3.24) 


1.50 x 10~ 3 (6.69 x 10~ 4 ) 


1.30 x 10" 1 (1.45 x 10" 2 ) 


10% 


16.78 (6.38) 


1.00 x 10" 3 (1.22 x 10" 3 ) 


2.30 x 10" 1 (3.42 x 10" 2 ) 



Table 1. Column(l): Noise level added to "observed" radial velocity map The quantities in Columns(2,3,4) were obtained by 
averaging over twenty one realisations for each level of noise. Column(2): pattern speed from using equation J7J, in units of 
kms _1 pc _1 . Columns(3,4): mean and rms residuals of recovered transverse velocities, respectively, in units of |fa;| max of the 
input map. 

and root-mean-squared (rms; o~r) value of the residual map were then calculated. When 
expressed in units of I^Lax of the input map, these were found equal to R = 1.09 x 1CT 3 
and o~r = 7.67 x 10~ 2 . These globally determined numbers should give the reader some idea 
of the dynamic range of the recovery method, when applied to noise-free spatially smoothed 
data. The spatial distribution of the errors in both the x and y velocities is best visualised 
with "arrow plots" of the velocity fields. Figure (2c) displays the velocity field of the input 
model in the disc plane. The reconstructed velocities are close to the model, and we do not 
present them separately Instead we plot the residual current field (recovered minus input) 
in Figure (2d). 

We also tested the recovery method on noisy data. In a real observation most of the 
error is likely to reside in the measurement of velocities, rather than the surface brightness. 
This is because the methods used to extract the velocity information from spectra are less 
robust than photometry. Therefore we added noise to Figure (lb), and kept Figure (la) 
noise-free. To each pixel of Figure (lb), we added Gaussian noise with mean equal to that 
observed, and rms equal to some fixed fraction of the mean. 3 . We experimented with three 
levels of noise, namely rms noise per pixel equal to 1%, 5%, and 10% of the mean. For each 
level of noise, twenty one realisations were explored. The pattern speed, x-current density 
and x-velocities were computed for each realisation, using equations (J7|) and (jSJ)- Comparing 
with the input model, we computed the mean (R) and rms (<tr) of the residual x-velocities 
for each realisation. The distributions of the twenty one i?'s and ctr's were peaked close to 
their mean values, R and cTr, respectively. The Table below provides estimates (and rms 
errors) for these, as well as for the pattern speed. 

The mean residual, R, is always quite small, implying that there is very little global sys- 
tematic shift in the x-velocities. This occurs because of cancellation between positive and 
negative residual velocities. The estimated pattern speed is also well-behaved, because this 

3 Adding noise to -F s ky is equivalent to adding noise to U = (-F s ky/2 s ky) , because we have kept £ s k y noise-free. 
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Figure 3. Recovered velocity fields and residuals for typical noisy realisations, (a) and (c) are the recovered velocity fields for 
noise levels of 5% and 10%, respectively, (b) and (d) are the respective residuals. 



is calculated using numbers from different Y cuts. However, the errors on Vl v increase dra- 
matically with noise, resulting in a significant increase in or . As earlier, the arrow plots are 
very revealing. The recovered and residual maps for the case of 1% noise are very close to 
the noise-free case discussed earlier. Therefore, in Figure 3 we display arrow plots only for 
noise levels of 5% and 10% . In addition to random errors in the residual velocities there 
are systematic alignments parallel to the line of nodes, the axis along which integrals are 
evaluated in the recovery method. However, as Figure 3 suggests, even for a noise level as 
high as 10%, the recovery method does not fail completely 
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4 CONCLUSIONS 



We have demonstrated that it is possible to recover the transverse velocities of steadily rotat- 
ing patterns in flat galaxies, using two dimensional maps of a tracer surface brightness and 
radial current density, if the tracer satisfies a source-free continuity equation. Our method is 
kinematical, and closely related to the TW method of determining pattern speeds. Indeed, 
the conditions that need to be satisfied — that the galaxy is flat, the pattern is steadily rotat- 
ing, and the tracer obeys a continuity equation — are identical to those assumed by TW84. 
Our main result is an explicit expression for the transverse velocities (equation |HJ), which is 
exact under ideal conditions. We have applied it successfully to simulated data, and demon- 
strated its utility in the presence of intrinsic numerical errors in the data, finite angular 
resolution, and noise. The TW relation for the pattern speed (equation EJ) emerges as an 
eigenvalue, and we expect our method to work well whenever the TW method gives a good 
estimate of the pattern speed. It is legitimate to be concerned that the conditions required 
to be satisfied might impose serious limitations in practice; the angle of inclination and line 
of nodes need to be estimated, the pattern need not be steadily rotating, the continuity 
equation need not be satisfied, the tracer distribution could be warped, the disc could be 
thick, and there could be streaming velocities in the z direction. All these are well-known 
worries about the applicability of the TW method itself. That they are not unduly restric- 
tive is evident from the suc cess that the TW method itself has enjoyed in the determination 
of pattern speeds (see e.g. 



1995 



Kent 



1987 



Gerssen. Kuiiken. and MerrifieL 



f in. J\u 
2001 :h 



kuiik en and T remaine 



1999; 



1991 



Merrifield and Kuiiken 



Bureau. Freeman. Pfitzner. and Meurei 



?urerlll99 
jlgOOa D( 



Baker 



2002; 



iDebattista and Willi 



Oers.se- 



Zimmer and Rand 



ams 



2001 



Debattis ta. Corsini. and Ascuerri 200 
Aguerri. Debattista. and CorsinilboO.^ . There 



Debattista. Gerha 



fore we are cautiously optimistic that our method of recovering transverse velocities can be 
applied usefully to two dimensional spectral maps. 
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